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1.  Introduction 

Thermoelectricity  is  actually  viewed  as  a  very  interesting  source 
of  electric  power  because  of  its  ability  to  convert  heat  flow  directly 
into  electricity.  In  particular,  thermoelectric  devices  as  energy  con¬ 
verters  are  easily  scalable  and  do  not  have  moving  parts,  or  liquid 
fuels.  These  excellent  features  make  them  applicable  in  almost  any 
situation  where  large  quantities  of  heat  tend  to  go  to  waste,  from 
clothing  to  large  industrial  facilities. 

Whilst  the  thermoelectric  materials  have  been  known  and 
understood  for  quite  some  time,  so  far  they  have  not  been  efficient 
enough  to  be  used  commercially.  Currently,  many  research  groups 
are  investing  their  efforts  in  finding  new  materials  with  appropri¬ 
ate  properties  to  better  use  thermoelectric  effects,  creating  so  effi¬ 
cient  thermoelectric  devices.  A  good  thermoelectric  material  will 
have  very  poor  thermal  conductivity  X,  but  a  very  high  electrical 
conductivity  ae,  so  that  the  so-called  figure-of-merit  Z  =  e2ae/X, 
with  e  being  the  Seebeck  coefficient,  which  is  the  main  responsible 
of  the  thermoelectric  efficiency,  is  as  high  as  possible. 

Indeed,  the  advent  of  nanotechnology  is  widening  the  range  of 
applicability  of  thermoelectric  materials.  In  fact,  the  nanostructures 
used  in  materials  maintain  good  electrical  conductivity,  while 
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reducing  the  thermal  conductivity.  The  performances  of  thermo¬ 
electric  devices  can  thus  be  enhanced  with  the  use  of  nanotechnol¬ 
ogy-based  materials  that  have  improved  thermoelectric  properties 
and  good  solar  energy  absorption  abilities.  Thermoelectric  materi¬ 
als  based  on  Bi2Te3,  for  example,  stand  out  as  perfect  examples 
highlighting  the  role  of  nanomaterials  for  thermoelectric  devices. 
Carbon  nanotubes  and  graphene  sheets  as  thermoelectric  materials 
also  exhibit  improved  thermoelectric  properties.  In  general,  the 
usage  of  nanostructures  smaller  than  the  wavelength  of  light 
enhances  the  scattering  of  photons  decreasing  so  the  thermal  con¬ 
ductivity.  This  decreasing  in  the  thermal  conductivity  seems  to  be 
the  most  vital  benefit  of  nanostructuring  for  thermoelectric 
materials. 

Although  it  is  very  clear  the  importance  of  using  nanotechnol¬ 
ogy  in  thermoelectricity,  the  design  of  good  thermoelectric  nan¬ 
odevices  is  still  far  from  its  optimal  solution.  This  principally 
because  the  physics  at  nanoscale  presents  some  dark  points,  as 
for  instance  the  role  played  by  memory,  nonlocal  and  nonlinear 
effects,  as  well  as  the  appropriate  definition  of  temperature  in  non¬ 
equilibrium  situations  [1-5], 

In  Ref.  [6]  we  have  investigated  the  influence  of  nonlocal  effects 
on  the  figure-of-merit  in  cylindrical  nanowires,  and  predicted  how 
it  depends  on  the  features  of  the  transversal  section.  In  the  present 
paper,  instead,  we  focus  our  attention  especially  on  the  role  of  the 
temperature  and  on  the  possibility  of  accounting  for  different  val¬ 
ues  for  electron  and  phonon  temperature.  Our  investigation  starts 
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Nomenclature 


cv  specific  heat  at  constant  volume 

tf  specific  heat  at  constant  volume  of  electrons 

cf  specific  heat  at  constant  volume  of  phonons 

E  electric-field  vector 

i  electric-current  density  vector 

J(s)  specific-entropy  flux  vector 

q  overall  local  heat-flux  vector 

q(c>  local  heat-flux  vector  due  to  electrons 

qlp)  local  heat-flux  vector  due  to  phonons 

s  specific  entropy 

T  average  temperature 

Te  electrons  temperature 

Tp  phonons  temperature 

u  specific  internal  energy 

ue  specific  internal  energy  of  electrons 

Up  specific  internal  energy  of  phonons 

Z  figure-of-merit 

Greek  symbols 

e  Seebeck  coefficient 

r]  thermoelectric  efficiency 

X  total  thermal  conductivity 

Xe  electrons  thermal  conductivity 


Xp  phonons  thermal  conductivity 

pe  electron  chemical  potential 

II  Peltier  coefficient 

p  mass  density 

g®  specific  electric  charge 

£  thermodynamic  state  space 

ae  electrical  conductivity 

c®  rate  of  entropy-density  production 

Te  relaxation  time  of  electrons 

T,-  relaxation  time  of  electric  charges 

xp  relaxation  time  of  phonons 

Subscripts 
e  electrons 

eff  effective 

i  electric-current  density 

max  maximum 

p  phonons 

Superscripts 
(e)  electrons 

(p)  phonons 

(s)  entropy 


by  the  observation  that,  as  the  electron  mean-free  path  4  is  usually 
shorter  than  the  phonon  mean-free  path  tv ,  when  heat  propagates 
in  a  system  whose  characteristic  size  L  is  such  that  4  <  L  <  (p,  it  is 
expected  a  very  high  number  of  electron  collisions,  and  only  scant 
phonon  collisions.  This  yields  that  the  electron  temperature  Te  may 
reach  its  local-equilibrium  value,  whereas  the  phonon  temperature 
Tp  is  still  far  from  its  own  local-equilibrium  value.  Conversely, 
when  the  electron  mean-free  path  (corresponding  to  the  elec¬ 
tron-phonon  collisions)  is  large,  one  may  have  the  so-called  phe¬ 
nomenon  of  “hot  electrons”,  namely,  a  population  of  electrons 
whose  average  kinetic  energy  (i.e.,  the  kinetic  temperature)  is  con¬ 
siderably  higher  than  that  of  the  phonons  [7,8],  Another  situation 
in  which  it  is  possible  to  have  different  temperatures  for  phonons 
and  electrons  is  when  a  high-frequency  electromagnetic  radiation 
is  used  to  supply  energy  to  the  system,  since  in  such  a  case  the 
electrons  receive  energy  at  a  rate  higher  than  that  at  which  they 
give  energy  to  the  phonons.  Therefore,  it  would  be  interesting  to 
find  possible  ways  to  measure  both  temperatures  [9,10], 

Besides  being  appealing  from  the  theoretical  point  of  view, 
accounting  for  two  different  temperatures  is  also  important  in 
practical  applications,  since  it  leads  to  a  more  realistic  computation 
of  the  thermodynamic  efficiency  of  the  thermoelectric  devices. 

Here  we  develop  a  mesoscopic  model  of  enhanced  thermoelec¬ 
tric  equations  which  account  for  different  phonon  and  electron 
temperatures.  Following  the  way  drawn  in  Refs.  [6,11],  we  also 
assume  that  the  overall  heat  flux  q  has  two  different  contributions: 
the  phonon  heat  flux  q®  and  the  electron  heat  flux  q(e),  in  such  a 
way  that  q  =  q(p)  +  q(f). 

The  layout  of  the  paper  is  the  following.  In  Section  2,  we 
develop  a  theoretical  model  describing  thermoelectric  effects 
when  the  electrons  and  the  phonons  do  not  have  the  same  temper¬ 
ature.  As  a  consequence  of  the  second  Kelvin  relation,  we  suggest  a 
possible  way  to  measure  those  temperatures.  In  Section  3,  as  a 
practical  application  of  that  model,  we  determine  the  efficiency 
of  a  thermoelectric  generator.  We  point  out  that  the  difference 
between  electron  and  phonon  temperature  can  contribute  to 
improve  the  thermoelectric  efficiency.  In  Section  4  we  draw  the 


main  conclusions  and  we  underline  that  our  model  may  also  cope 
with  the  phonon-drag  phenomenon,  suggesting  so  a  further  exper¬ 
iment  to  check  both  the  electron  and  the  phonon  temperature. 

2.  The  phenomenological  laws 

The  analysis  of  coupled  transport  processes  is  one  of  the 
outstanding  aspects  of  the  classical  theory  of  nonequilibrium  ther¬ 
modynamics  [3,12],  In  the  present  section  we  aim  to  derive  the 
phenomenological  laws  describing  the  thermoelectric  effects 
whenever  the  different  heat  carriers  (i.e.,  the  phonons  and  the  elec¬ 
trons  in  our  case)  no  longer  have  the  same  temperature. 

Since  we  assume  that  the  heat  carriers  behave  as  a  mixture  of 
gases  flowing  through  the  crystal  lattice  [13,14],  it  seems  logical 
to  suppose  that  the  internal  energy  of  phonons  per  unit  mass  up, 
the  internal  energy  of  electrons  per  unit  mass  ue  and  the  electrical 
charge  per  unit  mass  of  electrons  p(e)  belong  to  the  state  space.  In 
particular,  we  assume  that  those  state-space  variables  are  ruled  by 
the  following  evolution  equations: 

pdtUv  =  -V  •  qw  (la) 

pdtue  =  -  V  ■  q(e)  +  E  •  i  (lb) 

pdtQ(e)  =  -V  •  i  (lc) 

with  p  as  the  mass  density,  E  as  the  electric  field,  and  i  as  the  elec¬ 
tric-current  density  [6[.  We  note  that  the  summation  of  Eqs.  (la) 
and  (lb)  turns  out  the  well-known  energy-balance  equation 
pdtu  +  V  •  q  =  E  •  i 

obtained  in  Ref.  [15]  in  the  absence  of  a  magnetic  field,  once  the 
total  internal  energy  per  unit  mass  of  the  system  u  is  supposed  to 
be  given  by  the  constitutive  relation 

u  =  up  +  ue  (2) 

According  with  the  basic  principles  of  Extended  Irreversible 
Thermodynamics  (EIT)  [3,5],  we  may  assume  that  the  fluxes  of  pre¬ 
vious  unknown  variables  (namely,  q!pi,  q(e)  and  i)  are  the  other 
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state-space  variables.  Elevating  the  fluxes  to  the  status  of  indepen¬ 
dent  variables  amounts  to  introduce  memory  and  nonlocal  effects 
into  the  formalism  [16,17], 

In  the  Appendix  at  the  end  of  the  paper  it  is  shown  that,  when¬ 
ever  the  relaxation  times  of  those  fluxes  are  negligible,  the  choice 
above  of  the  state-space  variables  implies  that  the  specific  entropy 
s  is  such  that  s  =  s(up;  ue;  g(e>).  Therefore,  in  that  case,  from  the 
Gibbs  relation  between  s  and  the  state  variables,  we  have 


i  us  t  us  ,  US  , 

ds  =  —  dUp  +  —  due  +  dg(e) 
dup  p  due  dgW  u 

dtS  =  Y dtUp  +  j-dtUe  -  -gjjrfi 


wherein  Tp  =  ( ds/dup)~ 1 ,  Te  =  (ds/due)  1 ,  and  fie/ (g(e,Te)  =  -ds/dg(e) , 

jie  being  the  chemical  potential  of  the  electrons. 

Although  Eq.  (3)  does  not  present  any  problem  from  the  theo¬ 
retical  point  of  view,  from  the  practical  point  of  view,  instead,  it 
may  lead  to  some  perplexities  due  to  the  presence  of  Tp  and  Te 
therein.  In  fact,  one  may  naturally  wonder  whether  they  are  mea¬ 
surable  quantities,  or  not  [9,10],  Referring  the  readers  to  the  end  of 
the  present  section  for  a  first  possible  answer  to  that  question,  now 
let  us  only  observe  that  we  are  allowed  to  postulate  the  following 
further  constitutive  equations  which  relate  the  partial  internal 
energies  appearing  in  Eq.  (1 )  to  those  temperatures: 


Up  =  tfTp  (4a) 

Ue  =  C^Te  (4b) 

wherein  c(p)  and  c'p  are  the  phonon  and  the  electron  specific  heats 
at  constant  volume  [18],  respectively.  As  a  consequence  of  Eq.  (4), 
since  the  total  internal  energy  u  can  be  expressed  through  the  aver¬ 
age  temperature  T  as  u  =  c„T,  being  c„  =  c(p)  +  c(„e)  the  specific  heat 
at  constant  volume  of  the  whole  system  [19],  from  the  coupling  of 
Eqs.  (2)  and  (4)  we  obtain 


(5) 


which  states  a  very  strict  link  between  Tp,  Te  and  T,  the  latter  being 
a  measurable  quantity  in  practical  applications.  Note  that  in  the 
very  general  case  c'f'  and  c'f1  should  be  temperature-dependent 
functions,  but  here  we  deal  only  with  the  simplest  situation  in 
which  those  material  functions  are  constant,  in  order  to  emphasize 
the  essential  physical  ideas  and  their  consequences. 

In  few  words,  we  regard  the  phonons  and  electrons  as  a  mixture 
of  gases  flowing  through  the  crystal  lattice  [13,14],  each  of  which  is 
endowed  with  its  own  temperature.  In  this  way,  according  with 
the  theory  of  fluid  mixtures  with  different  temperatures  [20-22], 
we  are  allowed  to  assume  that  each  constituent  obeys  the  same 
balance  laws  as  a  single  fluid,  and  it  has  its  own  temperature. 
The  average  temperature  of  the  mixture  has  been  introduced  by 
the  consideration  that  the  internal  energy  of  the  mixture  is  the 
same  as  in  the  case  of  a  single-temperature  mixture  [21], 

The  substitution  of  Eq.  (1)  into  Eq.  (3)  leads  to 


V  •  q(p>  V  •  q(e) 


gWTe  V  ^  Te 
-_V  ,  q(C>  &  A  q(p)  ■  VTp  qW  ■  VTe 

\TP  +  T»  qWT  e)  Tp  T] 


Recalling  that  the  ti 
the  balance  law 


(i%)1VT'+T  (6) 

me  rate  of  the  specific  entropy  has  to  obey 


pdts  =  -V  •  J(s)  +  cr<s> 


(7) 


with  J(s)  being  the  entropy  flux,  and  a®  the  entropy  source,  the 
comparison  between  Eqs.  (6)  and  (7)  leads  to  the  following 
identifications: 


=  ^(“)-X(a)  (8b) 


wherein  J™  is  the  thermodynamic  flux,  and  X(IX)  is  its  conjugated 
thermodynamic  force  [3,12], 

Experience  indicates  that  Jw  and  X(a)  are  not  independent,  but 
that  there  exists  a  relationship  between  them.  Moreover,  it  has 
been  observed  that,  for  a  large  class  of  irreversible  processes,  the 
thermodynamic  fluxes  are  linear  functions  of  the  forces,  to  a  good 
approximation  [5,12,23],  This  observation,  which  is  also  the  sim¬ 
plest  way  to  ensure  that  a's)  is  a  non-negative  quantity  whatever 
the  thermodynamic  process  is,  allows  us  to  write  the  following 
phenomenological  relations  for  the  fluxes  appearing  in  our  theo¬ 
retical  model: 


-q^^2f+4-v(ife); 


w~q{e)  =l2'- 

i=ia^ 


(9a) 

(9b) 

(9c) 


In  Eq.  (9)  the  quantities  Lxp  mean  the  phenomenological  coeffi¬ 
cients,  which  are  related  to  experimental  quantities.  The  transport 
coefficients  are  related  to  the  thermo-physical  properties  of  the 
material  at  hand.  In  general,  the  depend  on  the  elements  of  the 
state  space,  as  well  as  on  the  thermodynamic  forces.  In  such  a  case 
we  are  dealing  with  a  nonlinear  nonequilibrium  theory  (NLNET) 
[24],  If,  instead,  these  coefficients  only  depend  on  the  elements 
of  the  state  space,  then  we  face  with  a  semi-linear  nonequilibrium 
theory  (SLNET)  [24],  Finally,  if  the  Lx!t  are  constant,  then  we  deal 
with  a  linear  nonequilibrium  theory  (LNET)  [24],  The  coefficients 
Laa  relate  the  thermodynamic  current  to  its  own  conjugated  ther¬ 
modynamic  force.  The  cross-coefficients  with  a.  #  /i,  instead, 
are  representative  of  the  coupling  between  different  physical 
effects  [12,25],  Any  phenomenon  in  which  two  or  more  transport 
effects  are  coupled,  such  as  thermal  and  electrical  conductivity, 
or  thermal  conductivity  and  diffusion,  is  called  cross-effect  [12,25], 
As  the  thermoelectric  effects  arise  from  the  physical  interrela¬ 
tion  between  heat  flow  and  electric  current,  in  the  classical  ther¬ 
moelectric  models  the  cross-coefficients  (with  ot  ¥=  /?)  are  only 
related  to  the  coupled  transport  of  heat  and  electricity.  In  the  pres¬ 
ent  model,  instead,  the  thermoelectric  effect  is  driven  by  three  gen¬ 
eralized  thermodynamic  forces,  namely,  the  phonon  and  electron 
temperature  gradient,  and  the  force  due  to  the  electric  field  and 
to  the  chemical  potential  of  the  electric  charge.  As  a  consequence 
the  cross-coefficients  I12  and  I2 1  account  for  the  cross  effects  due 
to  the  different  temperatures  of  the  heat  carriers  (which  are  lack¬ 
ing  in  the  classical  case),  the  coefficients  L,3  and  L3I  represent 
the  coupling  between  the  phonon  temperature  gradient  and  the 
electric  current,  and  the  coefficients  L23  and  L32  account  for  the 
coupling  between  the  electron  temperature  gradient  and  the  elec¬ 
tric  current.  In  the  standard  thermoelectric  models,  with  Tp  =  Te, 
we  have  only  two  cross-coefficients  and  a  2  x  2  transport  matrix. 
Here,  instead,  we  have  six  cross-coefficients  and  a  3  x  3  transport 
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matrix.  To  be  compatible  with  the  second  law  of  thermodynamics 
(expressed  as  the  positive-definite  character  of  the  entropy  produc¬ 
tion)  the  matrix  of  phenomenological  coefficients  must  be  posi¬ 
tive  definite  [12],  On  the  other  hand,  the  Onsager  reciprocal 
relations  [26,27]  ensure  that  the  transport  matrix  is  symmetric,  so 
that  the  Sylvester  criterion  concerning  the  positive  definiteness  of 
real-symmetric  n  x  n  matrices  is  applicable.  Such  a  criterion  states 
that  the  positiveness  of  all  the  leading  principal  minors  of  the  matrix 
is  necessary  and  sufficient  to  ensure  that  it  is  positive  definite  [28], 
Recalling  that  in  our  theory  q  =  qW  +  q(e>,  in  a  somewhat  differ¬ 
ent  form  the  phenomenological  equations  (9)  may  be  also  written  as 

q  =  -  {Ap  +  Aep )  VTP  -  {Ae  +  v)VTe  +  (jjfc  +  n )  i  (10a) 

§M-ffeeVre  +  <7e  |e  -  (10b) 

once  the  following  identifications  are  made: 

(  Ln  =  ApTp,  Li2  =  ApeTe ;  Lj3  =  0 

l  L21  =  AepTp:  L22  =  XeTe  +  <Te£Ten-  L23  =  -(I, IT  (11) 

(l3i=0;  L32  =  -aecTe:  L33  =  ae 

with  II  as  the  Peltier  coefficient,  and  Ap  and  Ae  as  the  thermal  conduc¬ 
tivities  of  the  material  whenever  the  sole  phonons  or  electrons  [6,29], 
respectively,  are  the  heat  carriers.  It  is  worth  observing  that,  as  the 
phonon  and  electron  internal  energies  (or,  due  to  Eq.  (4),  the  phonon 
and  electron  temperatures)  enter  the  state  space,  our  model  is  devel¬ 
oped  within  the  frame  of  a  semi-linear  nonequilibrium  theory. 

Indeed,  since  we  are  regarding  the  phonons  and  the  electrons  as 
constituents  of  a  same  mixture  of  flowing  heat  carriers,  then,  in 
principle,  one  has  to  take  into  account  the  possible  interactions 
between  them,  which  yield  further  thermal  contributions.  There¬ 
fore,  in  Eq.  (10a)  the  material  functions  Ape  and  7ep  express  those 
contributions  to  the  total  thermal  conductivity  of  the  material. 
Referring  the  readers  to  the  Appendix  for  more  comments,  here 
we  only  note  that  the  use  of  these  thermal  conductivities,  arising 
from  the  phenomenological  coefficients  L12  and  I2i  in  Eqs.  (9a) 
and  (9b),  respectively,  are  representative  of  the  cross  effects  in 
the  constitutive  equations  for  the  diffusive  fluxes  qW  and  q(f). 
However,  as  it  will  be  seen  in  Section  3,  these  cross  effects  do 
not  play  any  relevant  role  on  the  efficiency  in  the  thermoelectric 
energy  conversion,  being  the  difference  in  the  two  temperatures 
the  principal  responsible  for  possible  enhancements  of  it. 

In  Eq.  (11)  we  assumed  Li3  =L3i  =0  since  it  seems  logical 
because  the  phonons  are  not  expected  to  be  directly  sensitive  to 
the  external  electric  field,  at  least  in  a  first  approximation.  For  polar 
lattices  this  possibility  would  be  open,  and  Li3  could  be  different 
from  zero.  However,  here  we  take  the  simplest  expression. 

Owing  to  the  Onsager  reciprocal  relations  [26,27],  from  Eq.  (1 1 ) 
we  obtain 


Ape  -  Aep  (^) 
IT  =  eTe 


(12a) 

(12b) 


According  to  the  Sylvester’s  criterion,  the  additional  constraints 

Ap  >  0  (13a) 

ApiAe  +  ffeU)  >ApeAep  (13b) 

ApAe  >  ApeAep  (13c) 


are  necessary  and  sufficient  to  ensure  that  the  entropy  production  is 
non-negative  along  arbitrary  thermodynamic  processes.  In  fact,  due 
the  positive  definiteness  of  the  matrix  of  the  transport  coefficients, 
the  entropy  production  is  positive  whenever  at  least  one  thermody¬ 
namic  force  is  different  from  zero,  and  vanishes  when  all  the  ther¬ 
modynamic  forces  are  zero.  Such  a  situation  characterizes  the 


quasi-static  (i.e.,  reversible)  processes  [30],  which  correspond  to 
zero  entropy  production. 

The  relation  (12b)  is  well  known  in  thermoelectricity  when  Te  is 
replaced  by  T.  In  that  case  it  is  referred  to  as  the  second  Kelvin 
relation,  and  it  has  been  used  since  the  early  days  in  the  study  of 
thermoelectric  phenomena.  The  different  result  predicted  by  Eq. 
(12b)  with  respect  to  the  usual  statement  of  the  second  Kelvin  rela¬ 
tion  (i.e.,  II  =  eT)  should  not  be  considered  as  a  surprising  and 
unexpected  result.  In  fact,  the  classical  statement  is  correct  in  the 
case  that  phonons  and  electrons  have  the  same  temperature.  In 
the  case  they  have  different  temperatures,  instead,  the  relation 
(12b)  seems  more  appropriate  with  respect  to  the  classical  form. 
That  result  is  also  not  against  Onsager  reciprocity  relations 
[26,27]  but,  on  the  contrary,  it  is  an  illustration  of  them,  more  accu¬ 
rate  and  precise  than  the  classical  expression  with  a  single  temper¬ 
ature,  as  it  is  obtained  from  the  symmetry  of  coefficients  Lap 
appearing  in  Eq.  (11). 

Equation  (12b)  may  be  also  very  useful  in  practical  applications 
if  one  is  wondering  to  measure  Tp  and  Te.  In  fact,  by  means  of  a 
usual  thermometer,  one  is  only  able  to  check  the  average  temper¬ 
ature  T,  defined  by  Eq.  (5).  Equation  (12b)  turns  out  that  Te  should 
be  given  as  the  ratio  between  n  and  e.  Therefore,  if  one  is  able  to 
measure  n  and  e  then,  by  the  coupling  of  Eqs.  (5)  and  (12b),  in 
principle,  it  will  be  possible  to  determine  Tp,  too.  In  fact,  setting 

hcr  1-t^  n4) 

[h-b  k-Tf 

from  Eq.  (5)  we  have 

aft  +  (1  -  a) ft  =  1  «=*  ft  =T-Lr  (y^) ft  (15) 

Finally,  observing  that  a  6  ]0;  1[,  the  physical  constraint  ft  >  0 
implies  that  ft  e]0;a  1  [.  From  Eq.  (15)  it  is  easy  to  see  that  the 
condition  ft  —  1  =>  ft  —  1 ,  namely,  in  this  case  Te  =  TP  =  T:  the 
two-temperature  model  described  by  Eq.  (10)  reduces  to  the  usual 
single-temperature  model  [3,11],  i.e., 

i--,vr+(^+n)i 

i.-«vr  +  4-v(«L)] 

once  the  thermal  conductivity  A  of  the  material  is  supposed  to  be 
given  as  A  =  AP  +  Ae  +  2Ape.  Previous  considerations  allow  us  to 
claim  that  our  model  also  shows  a  strong  enough  physical  insight. 

3.  Efficiency  of  thermoelectric  generators 

Research  in  recent  years  has  been  focused  on  developing  both 
thermoelectric  structures,  and  materials  that  have  high  efficiency. 
In  the  present  section  we  point  out  the  influence  of  accounting  for 
two  different  temperatures  in  the  calculation  of  the  efficiency  r\  of 
a  thermoelectric  generator,  defined  as  the  ratio  between  the  elec¬ 
tric-power  output  Pei  and  the  total  heat  supplied  per  unit  time  Q, 
namely, 

<’6> 

To  this  end,  for  the  sake  of  simplicity,  we  consider  a  single 
one-dimensional  (y  is  the  sole  cartesian  coordinate)  thermoelectric 
element  of  length  L  under  steady  conditions.  The  hot  side  is  held  at 
a  temperature  Th  (assumed  to  be  the  upper  side,  at  y  =  I),  and  the 
cold  side  at  the  temperature  Tc  (the  lower  side,  at  y  =  0).  We  also 
assume  that  i  and  Q  enter  uniformly  into  the  hot  side  of  the  ele¬ 
ment.  In  such  a  situation,  the  electric-power  output  is 
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Pei  =  i  ■  Edy  =  ie(The  -lf)~ 


(17) 


once  the  electric  field  is  given  by  Eq.  (10b)  with  vanishing  values  of 
V  (yUe/g(e)),  for  the  sake  of  simplicity,  and  e  and  oe  do  not  depend  on 
the  temperature. 

From  Eq.  (10b),  instead,  we  obtain  that  the  total  heat  supplied 
per  unit  time  is 


4  =  Jg  qd y  = 


(jp  -  n)  (je  -  n) 

..A-L - V-La,  A„A - La 


(18) 


wherein  Ap  —  ).p  +  Xep,  and  A e  =  le  +  Zpe. 
Inserting  Eqs.  (17)  and  (18)  into  Eq.  (16) 


ie(The  -ifj-  i2Lo -1 

^  Ap  (Thp  -Tcpy-'  +  Ae  ( The  -lf)L-'  +  ni 


obtain 


(19) 


If  we  introduce  in  Eq.  (19)  the  coefficients  ft  and  ft  defined  in 
Eq.  (14),  and  take  into  account  Eq.  (12b),  by  straightforward  calcu¬ 
lations  we  have 


wherein  we  set  2  =  Ap  +  Ae,  y  =  (Ap/2)(ft/ft  -1).  and  x  is  the 
following  ratio  between  the  electric  current  and  the  heat  flux: 


■“HFCrJ  (2,) 

Moreover,  in  Eq.  (20)  r/c  =  1  -  Tc/Th  is  the  usual  Carnot 
efficiency,  and  i]r  is  a  reduced  efficiency.  Since  ijc  represents  the 
ideal  limit  of  the  thermodynamic  efficiency,  in  practical  applica¬ 
tions  one  should  find  the  right  way  to  enhance  rjr  in  order  to  have 
a  good  thermoelectric  efficiency.  Indeed,  it  is  easy  to  see  that 
whenever  the  ratio  x  defined  above  gets  the  value 

with  Z  =  e2(7e/2,  then  the  reduced  efficiency  gets  its  maximum 
value,  and  the  thermoelectric  efficiency  reads 


ft  <  1  =5-  Tp  >  Te,  ft  =  1  =>  Tp  =  Te,  and  ft  >  1  =>  Tr  <  Te. 


which  is  approximately  equal  to  the  electron  thermal  conductivity. 
Moreover,  we  supposed  that  ZT  =  1. 

As  it  can  be  seen,  the  maximum  efficiency  increases  for  increas¬ 
ing  values  of  ft.  This  means  that  the  bigger  Te  with  respect  to  Tp, 
the  better  the  performances  of  thermoelectric  devices.  Indeed, 
Fig.  1  also  allows  to  analyze  the  role  played  by  c'f1  and  cf\  the  lat¬ 
ter  being  usually  higher  than  the  former.  In  fact,  it  points  out  that 
the  (positive)  slope  of  the  curve  tymax  =  z/max(j8 1)  gradually 
decreases  whenever  a  assumes  very  small  values,  whereas  it  basi¬ 
cally  takes  a  constant  value  when  a  reaches  high  enough  values,  in 
such  a  way  that  whenever  Te  >  Tp,  the  higher  a,  the  higher  f?max. 

At  the  very  end,  we  observe  that  in  Fig.  1  the  value 
>1mlx/*lc  -  0.17,  attained  whenever  ft  =  1 1  corresponds  to  the  case 
of  a  single-temperature  model.  Thus  we  may  conclude  that 
whenever  the  electron  temperature  is  higher  than  the  phonon 
temperature,  our  two-temperature  model  yields  an  efficiency  which 
is  higher  than  that  of  the  usual  single-temperature  model. 

4.  Conclusions 


ZTft+  2(y  +  l)(^l-/l+ffJ 
Wx 


(23) 


which  reduces  to  the  classical  form  for  the  maximum  thermoelec¬ 
tric  efficiency  [3,31]  whenever  Tp  and  Tr  coincide,  i.e.,  when 

ft  =  ft  =  1  and  y  =  0. 

From  Eq.  (23)  it  easy  to  recover  the  usual  result  that  the  larger 
the  figure-of-merit,  the  higher  the  efficiency  of  a  thermoelectric 
device. 

However,  Eq.  (23)  clearly  points  out  that  also  the  differences 
between  Tp  and  Tc  influence  rjmax,  whereas  the  cross  effects  related 
to  the  phenomenological  coefficients  I]2  and  I2i  in  Eqs.  (9a)  and  (9b) 
substantially  remain  unaltered  the  predictions  of  our  two-tempera¬ 
ture  model  with  respect  to  those  of  the  usual  single-temperature 
model.  This  may  be  interesting  in  practical  applications,  since  the 
most  part  of  the  research  groups  is  focusing  the  attention  only  on 
the  search  of  new  materials  with  high  values  of  the  figure-of-merit. 

In  Fig.  1  we  plot  the  behavior  of  the  ratio  »7max/t7c  as  a  function  of 
ft  for  two  different  values  of  the  nondimensional  parameter  a,  i.e., 
a  =  0.05  and  a  =  0.5.  For  the  sake  of  illustration  in  our  computa¬ 
tion  we  assumed  Ap  =  Ae,  as  the  materials  commonly  used  in 
thermoelectric  applications  show  a  phonon  thermal  conductivity 


Thermoelectric  devices  have  gained  importance  in  recent  years 
as  viable  solutions  for  appealing  applications  such  as  spot  cooling 
of  electronic  components,  remote  power  generation  in  space  sta¬ 
tions  and  satellites,  etc.  These  solid-state  devices  have  long  been 
known  for  their  reliability,  rather  than  their  efficiency.  They  do 
not  contain  moving  parts,  and  their  performances  primarily  rely 
on  the  material  selection  which,  indeed,  has  not  yet  achieved 
sufficiently  satisfactory  candidates  for  widespread  practical  appli¬ 
cations.  Therefore,  the  research  in  recent  years  has  been  focused 
on  developing  both  thermoelectric  structures,  and  materials  that 
have  higher  efficiencies  than  that  of  the  currently  available  mate¬ 
rials.  The  advent  of  nanotechnologies,  which  allow  to  improve 
the  actual  performances  of  thermoelectric  devices,  has  shifted 
the  focus  towards  a  deeper  understanding  of  the  carrier-transport 
behavior  in  nanostructures.  In  fact,  the  quantum  confinement  in 
nanostructures  increases  the  local-carrier  density  of  states  per  unit 
volume  near  the  Fermi  energy,  yielding  an  enhancement  in  the 
Seebeck  coefficient  [32],  while  the  thermal  conductivity  can  be 
decreased  due  to  phonon  confinement  and  phonon  scattering 
[33,34],  This  way,  the  combined  benefits  of  reduced  thermal 


1  Along  with  previous  observations,  in  this  case  one  also  has  ft  =  1.  and  moreover 
7  =  0. 
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conductivity  and  improved  Seebeck  coefficient  imply  theoretically 
higher  values  of  ZT  compared  to  the  bulk  structures.  However, 
experimental  observations  have  not  been  able  to  achieve  the  pre¬ 
sumed  benefits  of  nanostructured  thermoelectric  devices,  despite 
theoretically  predicted  improvements. 

Owing  to  the  need  for  a  better  understanding  of  the  effect  of  all 
the  significant  factors  contributing  to  the  thermoelectric  figure-of- 
merit  of  nanoscale  devices,  in  a  previous  paper  [6]  we  have 
investigated  the  influence  of  nonlocal  effects  on  it  in  cylindrical 
nanowires,  and  we  predicted  how  Z  depends  on  the  features  of 
the  transversal  section. 

Here,  instead,  we  have  analyzed  the  possible  consequences  of 
accounting  for  different  temperatures  for  the  different  heat  carriers 
which  may  be  interesting  in  practical  applications,  as  we  pointed 
out  in  Sections  1  and  3.  In  particular,  in  the  framework  of  EIT 
[3,5],  in  the  present  paper  we  have  developed  a  model  of  enhanced 
thermoelectric  equations  assuming  that  the  population  of  heat  car¬ 
riers  behaves  as  a  mixture  of  flowing  gas  particles  [6,13,14]  with 
different  temperatures,  according  with  the  theory  proposed  by 
Ruggeri  and  coworkers  [20-22],  It  seems  worth  noticing  that 
accounting  for  different  temperatures  is  also  a  well-known  mech¬ 
anism  of  heat  conduction  in  complex  materials,  where  one  of  the 
components  can  be  thermally  excited  independently  of  the  other 
[7,8,35,36],  In  the  case  of  thermoelectric  effects,  instead,  we  feel 
that  the  consideration  of  two  temperatures  is,  at  least,  a  theoretical 
valuable  approach  because  it  allows  to  identify  with  better  concep¬ 
tual  clarity  which  part  of  the  several  effects  is  directly  related  to 
electrons  (or  holes),  which  one  to  phonons,  and  which  one  to  their 
mutual  interactions.  Thus,  it  is  worth  to  explore  it  for  its  own  con¬ 
ceptual  sake,  as  a  model  wider  than  the  usual  one,  working  directly 
with  a  single  temperature.  On  the  other  side,  this  may  open  the 
possibility  of  new  practical  strategies,  although  the  effects  are 
expected  to  be  small  in  usual  circumstances.  In  fact,  the  results 
in  Eq.  (23),  arising  from  the  two-temperature  model  in  Eq.  (10), 
point  out  that  the  bigger  the  electron  temperature  than  the  phonon 
temperature,  the  higher  the  performances  of  a  thermoelectric 
device.  This  can  be  clearly  seen  in  Fig.  1  wherein  t/max  increases 
for  enhancing  values  of  the  nondimensional  ratio  ft  =  Te/T.  Figure 
1  also  allows  us  to  make  a  comparison  between  the  theoretical  pre¬ 
dictions  of  a  single-temperature  model  (ft  =  1 )  and  those  of  our 
two-temperature  model  (ft  ¥=  1 ).  In  particular,  whenever  ft  is 
higher  than  Tp,  it  points  out  that  the  latter  model  yields  a  larger 
efficiency  in  the  thermoelectric  energy  conversion  with  respect 
to  that  predicted  by  the  former  model. 

Since  from  the  theoretical  point  of  view  it  seems  very  clear  the 
importance  of  accounting  for  ft  and  ft,  it  would  be  very  useful  in 
practical  applications  to  find  a  possible  way  to  estimate  those  tem¬ 
peratures  [9,10],  To  this  end,  by  observing  that  our  theoretical 
model  predicts  that  the  classical  second  Kelvin  relation  breaks 
down,  as  stated  by  Eq.  (12b),  we  also  proposed  a  possible  way  of 
checldng  both  temperatures.  In  more  details,  we  observed  that  if 
one  is  able  to  measure  at  nanoscale  the  Peltier  and  the  Seebeck 
coefficients  separately,  then  their  ratio  would  turn  out,  in  principle, 
the  electron  temperature  and  then,  as  a  consequence  of  Eq.  (5),  the 
phonon  temperature,  too.  It  is  worth  to  note  that  at  nanoscale  the 
different  material  functions  may  deviate  from  their  corresponding 
bulk  values,  in  such  a  way  that  several  methods  of  measuring  them 
may  be  found  in  literature  [37,38], 

From  the  practical  point  of  view,  the  way  of  having  a  truly  dif¬ 
ferent  temperature  for  electrons  (and  holes)  and  lattice  (phonons) 
is  by  means  of  a  pulse-laser  excitation  yielding  its  energy  to  the 
charged  particles  (electrons),  which  later  share  their  energy  excess 
with  the  lattice.  Despite  the  lattice  may  be  also  charged,  the  much 
higher  mass  of  the  ions  makes  that  the  electrons  may  absorb  more 
energy  from  the  pulse.  However,  this  strategy  is  not  of  interest  for 
practical  thermoelectric  devices,  because  it  would  require 


spending  energy  on  the  laser,  from  which  only  a  part  would  be 
taken  by  the  system.  A  different  strategy  would  be  by  using  a  sys¬ 
tem  composed  of  two  (or  several)  thin  layers  (as  for  instance  a 
superlattice)  which  are  shorter  than  the  phonon  mean-free  path, 
but  larger  than  the  electron  mean-free  path.  If  the  electron  and 
phonon  contributions  to  the  heat  flux  within  the  layers  are,  respec¬ 
tively,  q(e)  =  -ftVft  and  q'pi  =  -ftVft,  and  the  temperature  dis¬ 
continuities  at  an  interface  are  Aft  =  Reqe  and  A Tp  =  Rpqp,  Re  and 
Rp  being  the  respective  thermal  resistance  of  the  interface,  one 
may  obtain  the  profile  for  ft  and  Tp  along  the  system.  Both  temper¬ 
atures  must  be  equal  to  the  heat  baths  at  the  two  ends  of  the  whole 
system,  but  they  may  be  different  from  each  other  along  the  sys¬ 
tem.  By  using  suitable  materials  to  have  ft.,  ft,  ft  and  ft  (or  a  suit¬ 
ably  reduced  effective  phonon  thermal  conductivity),  one  could 
have  regions  with  ft  higher  than  Tp,  and  regions  with  ft  lower  than 
ft.  A  detailed  analysis  of  both  profiles  should  be  carried  out  to 
study  the  global  effects  of  these  differences  on  the  efficiency  of 
the  thermoelectric  conversion  in  the  whole  device.  Of  course,  elec¬ 
trons  exchange  energy  with  phonons  everywhere.  Thus,  it  would 
be  convenient  that  the  electron  heat  flux  is  relatively  large,  in  order 
that  this  energy  exchange  does  not  bring  to  zero  the  temperature 
differences.  These  effects  should  be  studied  in  detail.  Though  the 
final  practical  outcome  may  turn  to  be  small,  the  new  detailed 
physical  understanding  of  the  system  would  be  worth  of  the  effort. 

Indeed,  in  our  two-temperature  model  the  problem  of  measur¬ 
ing  ft  and  ft  may  be  also  related  to  the  phonon-drag  phenomenon 
[39,40],  The  classical  theory  of  thermoelectricity,  in  fact,  is  based 
on  the  assumption  that  the  flow  of  charge  carriers  and  phonons 
can  be  treated  independently.  Under  this  assumption,  the  Seebeck 
coefficient  depends  solely  by  the  spontaneous  electron  diffusion. 
However,  when  the  two  flows  are  linked,  the  effect  of  electron- 
phonon  scattering  should  be  taken  into  account.  Hence,  in  general, 
the  Seebeck  coefficient  shows  two  independent  contributions:  the 
conventional  electron-diffusion  contribution  and  the  phonon-drag 
contribution  [40],  The  diffusion  part  is  caused  by  the  spatial  varia¬ 
tion  of  the  electronic  occupation  in  the  presence  of  a  thermal 
gradient,  whereas  the  drag  part  arises  by  the  interaction  between 
anisotropic  lattice  vibrations  and  mobile  charge  carriers.  The 
overall  phonon-drag  effect  leads  to  an  increase  in  the  Seebeck  coef¬ 
ficient.  If  we  look  at  Eq.  (10b),  we  may  observe  that  it  allows  to 
introduce  the  following  effective  Seebeck  coefficient 

eefr-eft^ft=M  (24) 

which  also  allows  to  claim  that  the  deviation  of  the  effective 
Seebeck  coefficient  from  its  bulk  value  represents  a  further  possible 
measurement  of  the  electron  temperature.  Moreover,  along  with 
previous  observations  about  the  phonon-drag  phenomenon,  from 
Eq.  (24)  we  further  claim  that  in  general  ft  should  be  greater  than 
unit,  namely,  ft  >  ft. 
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Appendix  A.  Compatibility  with  second  law  of  thermodynamics 


The  results  we  derived  in  the  present  paper  are  based  upon  the 
system  of  Eq.  (10).  In  Section  2  we  already  pointed  out  a  sufficient 
set  of  thermodynamic  conditions  ensuring  the  compatibility  of 
those  equations  with  the  second  law.  Indeed,  Eq.  (10)  follow  from 
the  assumption  that  the  specific  entropy  only  depends  on  the  basic 
variables  up,  ue  and  g,e\  whereas,  according  with  the  general 
principles  of  EIT  [3,5,16,17],  we  have  assumed  that  both  those  vari¬ 
ables,  and  the  independent  fluxes  q(p),  q(e)  and  i  belong  to  the  state 
space.  Here  we  spend  more  comments  about  that  assumption.  In 
doing  this,  we  take  advantage  of  the  second  law  of 
thermodynamics. 

The  second  law,  accounting  for  the  natural  evolution  of  a 
system  in  any  possible  thermodynamic  process,  in  fact,  is  a  very 
valuable  tool  to  study  the  physical  consistency  of  any  thermody¬ 
namic  model.  In  terms  of  the  entropy  density  s,  its  local  form  reads 
as  in  Eq.  (7),  wherein  the  entropy  source  has  to  be  non-negative  in 
any  admissible  thermodynamic  process.  Since  in  Eq.  (7)  neither  s, 
nor  J,s)  belong  to  the  state  space,  then  they  have  to  be  expressed 
through  constitutive  equations.  If  we  assume  for  them  the  very 
general  forms 


S  =  s(up:Ue;e<eW>p);q(e);i) 

J(s)  =  J(s)  (up  ;  ue ; !? (e) ;  q<p> ;  q(e) ; ') 


(25) 


then,  by  the  chain  rule,  from  Eq.  (7)  we  get 


wherein  tp,  te  and  t,  are  the  relaxation  times  of  phonons,  electrons 
and  electric  current,  respectively  [6,42],  In  steady-state  situations, 
or  whenever  those  relaxation  times  are  vanishing,  from  Eq.  (27)  it 
is  easy  to  recover  Eq.  (10). 

Before  to  go  further  in  the  thermodynamic  analysis,  we  feel  that 
Eq.  (27)  deserve  some  comments.  They  introduce  a  theoretical 
model  for  thermoelectric  effects  which  allows  to  point  out  clearly 
which  part  of  the  several  effects  is  directly  related  to  the  different 
heat/electric  carriers,  and  which  one  to  their  mutual  interactions. 
Mutual  interactions,  in  particular,  in  Eqs.  (27a)  and  (27b)  are  intro¬ 
duced  by  cross  effects  of  the  form  ApeVTe  and  AepVTp,  arising  from 
the  phenomenological  coefficients  I12  and  I2i  in  Eq.  (9)  which  are 
peculiar  of  our  two-temperatures  model.  We  already  observed  that 
these  effects  do  not  play  any  direct  relevant  role  on  the  thermo¬ 
electric  efficiency.  They  are  very  important,  instead,  if  one  is  won¬ 
dering  to  determine  both  the  Tp  profile  and  the  Te  profile  in  our 
system.  In  fact,  since  they  couple  the  differential  equations  for 
qw  and  q(e),  then  also  small  perturbations  in  one  temperature  have 
repercussions  on  the  other  temperature.  However,  we  note  that  the 
analysis  of  the  possible  influence  of  cross  effects  in  the  constitutive 
equations  for  thermoelectricity  is  a  very  interesting  research  play¬ 
ground,  since  it  may  allow  to  discover  new  ways  to  enhance  the 
performances  of  thermoelectric  devices,  as  it  has  been  pointed 
out  in  Ref.  [11],  for  example. 

Equations  (27)  follow  from  the  observation  that  each  thermo¬ 
dynamic  flux  Jla)  is  described  by  a  generalized  transport  equation 
of  the  form  [5] 

T*dJ(a)  +  f>  =  Up  ■  X(« 

wherein  Up  is  the  matrix  of  the  phenomenological  coefficients.  It  is 
worth  observing  that  whenever  the  relaxation  times  tp,  te  and  t, 
are  negligible,  the  summation  of  Eqs.  (27a)  and  (27b)  turns  out 
Eq.  (10a),  whereas  Eq.  (27c)  reduces  to  Eq.  (10b). 

To  derive  the  thermodynamic  restrictions  imposed  by  the 
inequality  >  0,  it  is  possible  to  follow  the  procedure  proposed 
by  Liu  [43]  which  introduces  the  constraints  in  Eqs.  (1 )  and  (27)  by 
means  of  suitable  Lagrange  multipliers,  in  such  a  way  that  the  fol¬ 
lowing  inequality  has  to  be  fulfilled: 


6q® 


4-jfiL 

dq<P> 


Vq® 


y:Vq'e> 


:  Vi 


(26) 


wherein  the  colon  stands  for  the  double  inner  product  of  two  matri¬ 
ces.  Second  law  of  thermodynamics  dictates  that  the  right-hand 
side  of  Eq.  (26)  must  be  nonnegative  for  arbitrary  thermodynamic 
processes.  It  contains  the  higher  derivatives  dtup,  dtue,  dtg(e>, 
3tq(e),  dti,  Vup,VuE,  Vqli,),  Vq(e)  and  Vi,  which  can 

assume  completely  arbitrary  values  [41],  as  they  account  for  the 
natural  and  independent  evolution  in  the  space-time  of  the  basic 
variables  up,  ue  and  Q[e\  as  well  as  for  the  evolution  of  the  indepen¬ 
dent  fluxes  q®,  q(e>  and  i. 

The  evolution  equations  of  up,  ue  and  o<e)  are  given  by  Eq.  (1), 


a*} 

(28) 


(7<s)  -  yp  ( pd,uv + V  •  qW)  -  y,  ( pd,ue  +  V  •  qw  -  E  •  i) 

-  y,  (pdtQ&  +  V  ■  i)  -  rp  (TjAqC”  +  q<”>  +  APVTP  +  V VTe) 

-II  |  Ted,q(e>  +  q(£>  +  V VTP  +  (Je  +  <JeeU) VTe  -  aeU  |e  -  V  ^ 

-r({rIa(i+i+ffcevrc-ae[E-v(^)]}>o 

Since  the  Lagrange  multipliers  yp,  ye,  yt,  rp,  re  and  r,  are 
well-defined  on  the  state  space  [41  ],  as  a  consequence  of  the  cou¬ 
pling  of  Eq.  (26)  and  inequality  (28)  thermodynamic  restrictions 
ensue  by  nullifying  the  coefficients  of  the  higher  derivatives  [41], 
In  particular,  setting  equal  to  zero  each  coefficient  of  the  time- 
derivatives  of  the  state-space  variables,  one  has: 


ds 


1 


—  =  yp  *=>  r  =  yp 


whereas  the 

time  rates  of  q®,  q(e)  and  i  in 

our  model  read, 

ds 

He 

respectively, 

dQW  ~  7i  ^ 

gWTe  7i 

Ti,atq(p)  +  q(p) 

m-ApVTp  -  APeVTe 

(27a) 

ds  rpxp 

Tedtq(e)  +  q(e) 

m^AepVTp  -  (Ae  +  <7een)vre 

dqM  p 

ds  Texe 

+  aenE-vfe)|+(j)i 

(27b) 

9q<e)  ~  P 

L  Vel,/J  WV 

ds  _  r,T,- 

Ttdti  +  i  =  -aeeVTe  +  ae  |e  -  V 

(27c) 

di~  p 

(29a) 

(29b) 

(29c) 

(29d) 

(29e) 

(29f) 
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These  relations  turn  out  useful  information  about  the 
dependence  of  the  specific  entropy  on  the  state-space  variables. 
A  possible  form  of  s  compatible  both  with  the  principle  of 
maximum  entropy  at  the  equilibrium,  and  with  Eq.  (29)  is 


s=s0(up-ue-Q^)  -- 


2{K  +  en  (je)T2e 


- ^ii 

2eoeT] 

provided  the  following  identifications  are  made: 


p  n[e) 

'  2(Ae  +  eUae)T2e  q 

i 

2eaeT2e 


'  q(e) 


(30) 


(31a) 

(31b) 

(31c) 


From  Eq.  (30)  it  follows  that  whenever  the  relaxation  times 
tp,  Te  and  t,  are  negligible,  then  the  dependence  of  s  on  the  fluxes 
qw,  q<e>  and  i  may  be  neglected  as  well,  and  one  may  assume 
s  rs  So  (Up ;  ue;  £(e)).  The  result  above  proves  that  the  Gibbs  equation 
(3),  postulated  in  Section  2,  is  in  accordance  with  second  law  of 
thermodynamics. 

Moreover,  if  one  sets  equal  to  zero  the  coefficients  of  the  first- 
order  spatial  derivatives  of  the  state-space  variables,  the  following 
further  set  of  thermodynamic  restrictions  ensues: 


I»  =  Wtr,[i,«,nA(i)]  +  IW<s;(&)  <32a> 
— -  L',  +  +  (  ((;•  )]} 


+r'4+a 

m 

(32b) 

£=".(nr.  +  r,)5|5| 

m 

(32c) 

=  y  U 

flqW  yp 

(32d) 

(32e) 

(321) 

with  U  being  the  unitary  matrix. 

It  is  easy  to  verify  by  direct  substitution  that  the  system  of 
thermodynamic  restrictions  (32)  admits  as  solution  the  constitu¬ 
tive  Equation  (8a). 

These  thermodynamic  considerations  allow  to  claim  that  the 
model  in  Eq.  (10),  which  are  a  particular  case  of  the  more  general 
Eq.  (27),  are  in  accordance  with  the  basic  principles  of  continuum 
physics. 
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